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We study the entropy of chiral 2+1 -dimensional topological phases, where there are both gapped bulk ex- 
citations and gapless edge modes. We show how the entanglement entropy of both types of excitations can 
be encoded in a single partition function. This partition function is holographic because it can be expressed 
entirely in terms of the conformal field theory describing the edge modes. We give a general expression for 
the holographic partition function, and discuss several examples in depth, including abelian and non-abelian 
fractional quantum Hall states, and p + ip superconductors. We extend these results to include a point contact 
allowing tunneling between two points on the edge, which causes thermodynamic entropy associated with the 
point contact to be lost with decreasing temperature. Such a perturbation effectively breaks the system in two, 
and we can identify the thermodynamic entropy loss with the loss of the edge entanglement entropy. From 
these results, we obtain a simple interpretation of the non-integer 'ground state degeneracy' which is obtained 
in 1 + 1-dimensional quantum impurity problems: its logarithm is a 2+1 -dimensional topological entanglement 
entropy. 



I. INTRODUCTION 

Entanglement is one of the characteristic peculiar features 
of quantum mechanics. It lay at the heart of the debates be- 
tween Bohr and Einstein, Podolski, and Rosen, and it is es- 
sential for the Bell's inequality violation which distinguishes 
quantum mechanics from classical hidden variables theories. 
Entanglement between logical qubits is a resource for quan- 
tum computation; entanglement between physical qubits and 
the environment is a threat to quantum computation; entangle- 
ment between physical qubits is essential for error correction. 

Entanglement entropy is one measure of entanglement. If 
a system can be subdivided into two subsystems A and B, 
then even if the whole system is in a pure quantum state |\&), 
subsystem A will be in a mixed state with density matrix p 
obtained by tracing out the subsystem B degrees of freedom: 

^=tr B (|*)<*|) (1) 

The entropy of this density matrix will be zero if the state | *f?) 
is the direct product of a pure state for subsystem A with a 
pure state for subsystem B or, equivalently, if the matrix p A 
has only a single non-zero eigenvalue. If the entropy 

S A = -ti(p A lnp A ) (2) 

is non-zero, then the subsystems A and B are entangled. 

The entanglement entropy Sa is a measure of the correla- 
tions between the degrees of freedom of the A and B subsys- 
tems. In a 1+1-dimensional quantum system, the dependence 
of Sa on the length L of subsystem A can be used to dis- 
tinguish critical from non-critical systems: it remains finite 
as the length increases, except at criticality, where it diverges 
logarithmically with a universal coefficient^. Of course, in 
this case, there are other measures of criticality, such as the 
power-law decay of correlation functions. In a gapped sys- 
tem, the leading term in the entanglement entropy is propor- 
tional to the surface area of the boundary. The coefficient is 



cutoff-dependent. However, in a 2 + 1-dimensional system in 
a gapped topological phase, the first subleading term is uni- 
versal and independent of the size or shape of A*&: 

S A = aL - InV + ... (3) 

It is not immediately obvious that T> in the above equation 
is uniquely defined since a is a cutoff-dependent coefficient 
(in two spacetime dimensions, the 'surface area' is the length 
of the boundary). However, by dividing a system into three or 
more subsystems and forming an appropriate linear combina- 
tion of the resulting entanglement entropies, the length term 
can be canceled, leaving only the universal term 3 4 . Such a 
construction can be used to give a more precise definition of 
the topological entanglement entropy, but its essential mean- 
ing is captured by 0. The quantity T> is the total quantum 
dimension of the topological phase, which we define in sec- 

tionini 

This is potentially a very useful probe of a topologi- 
cal phase. In such a phase, all correlation functions are 
topologically-invariant at distances longer than some finite 
correlation length and energies lower than a corresponding 
energy scale. Hence, in order to identify such a state, one 
must examine the ground state degeneracy on higher-genus 
surfaces or the braiding properties of quasiparticle excitations. 
The entanglement entropy gives us a handle on this topolog- 
ical structure stemming from the ground state wavefunction 
alone. 

Topological quantum field theories (TQFTs) were widely 
studied in the string-theory community in the late '80s and 
early '90s, following Witten's celebrated work 5 on Chern- 
Simons field theory, conformal field theory (CFT), and the 
Jones polynomial of knot theory 6 . There has been a recent 
revival of interest in TQFTs due to the possible emergence 
of topological phases in electronic condensed-matter systems 
and the proposed use of such phases as a platform for fault- 
tolerant quantum computation 7 ^. The above definition of a 
topological phase can be compactly (though tautologically) 
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restated as a phase of matter for which the long-distance ef- 
fective field theory is a topological field theory. The entangle- 
ment entropy is a way of extracting one of the basic parame- 
ters of a topological field theory, its quantum dimension. 

A chiral TQFT in 2+1 dimensions is related to 1+1- 
dimensional rational conformal field theory (RCFT) in several 
ways. When the TQFT is defined on a manifold with bound- 
ary, the boundary degrees of freedom form an RCFT. The 
topological invariance of the theory means that there is no dis- 
tinction between spacelike and timelike boundaries. Hence, 
the RCFT describes both (a) the ground state wavefunction(s) 
at some (2+0)-d equal-time slice when spacetime is M x /, for 
a compact surface M and time interval /, and (b) the dynam- 
ics of the (l+l)-d boundary when spacetime is X x R, where 
X is a surface with boundary. In the RCFT, the topological 
braiding/fusion structure of the TQFT is promoted to a holo- 
morphic structure with branch cuts implementing non-trivial 
monodromies. This structure occurs in the TQFT ground 
state(s) (the (2+0)-d case) as a result of the choice of holo- 
morphic gauge. In the (l+l)-d case, however, it represents 
the actual critical dynamics of the excitations of the edge of 
the system. 

Although other possible settings are also suspected, there is 
only one place in nature where topological phases are known 
for certain to exist: the quantum Hall regime. In this regime, 
this 'holographic' relationship between the (2+l)-d TQFT de- 
scribing the bulk and the (l+l)-d RCFT describing edge ex- 
citations is beautifully realized. Since it is easier to experi- 
mentally probe the edge, the topological properties of quan- 
tum Hall states have, thus far, been mostly probed through the 
constraints that they place on the dynamics of the edge. 

One can continue to still lower dimensions and consider 
boundaries or defect lines in 2d classical critical systems and 
impurities in (l+l)-d quantum critical systems with dynami- 
cal critical exponent z = 1. (In many (2+l)-d or (3+l)-d crit- 
ical systems, only the s-wave channel interacts with the impu- 
rity. By performing a partial wave decomposition and keeping 
only this channel, one can map the problem to a (l+l)-d prob- 
lem. Hence, the restriction to (l+l)-d is not very severe.) In 
these situations, the constraints of conformal invariance in the 
bulk strongly constrain the low-energy behavior of the bound- 
ary/defect correlations or impurity dynamics. The same meth- 
ods can be applied to a point contact in a (2+l)-d system in 
a topological phase. The point contact allows tunneling be- 
tween the gapless excitations at the edge and can, therefore, 
be understood as an impurity in a (l+l)-d system. 

In particular, we can write the entropy of the impurity prob- 
lem as 

5 imp = f L + Hg) + ■■■ (4) 

where L is the size of one-dimensional space. The leading 
piece is non-universal, but a universal part of the entangle- 
ment entropy proportional to the central charge of the confor- 
mal field theory can be extracted 9 . In this paper, we will focus 
on subleading pieces of the entropy like ln(g). Although ln(g) 
is often called the "boundary entropy", we see that this inter- 
pretation is not always useful: we will see how such a term can 



be present even when there is no boundary to the 1+1 dimen- 
sional system. Moreover, if it were truly the entropy of the 
impurity, then g would be an integer. However, if one com- 
putes it by first taking the thermodynamic limit LT — > oo, and 
then taking T — ► 0, g is not always an integeriSiiLi^ii^. This 
makes this subleading contribution to the entropy somewhat 
mysterious to interpret physically. 

To make this clearer, we will examine in detail the situation 
in which the (l+l)-d critical system is the edge of a (2+1)- 
d system in a chiral topological phase. In particular, we will 
show that the value of g in @ is closely related to that of 
the quantum dimension T> of the topological quantum phase, 
that occurs in the entanglement entropy in 0. Although they 
are both subleading terms in entropies, their relation is far 
from obvious. Indeed, their very definitions are quite differ- 
ent, involving an interchanged order of limits. The entangle- 
ment entropy is a property of the ground state obtained in the 
zero temperature limit, when the region size L is subsequently 
taken to infinity. On the other hand, extracting the subleading 
ln(g) correction of the thermodynamic entropy as in @ re- 
quires taking L to infinity before taking the zero temperature 
limit. Despite this opposite ordering of limits, we will show 
that both quantities follow from the same deep results in con- 
formal field theoryi2iiiiI£. 

The connection can be understood heuristically by consid- 
ering the entropy loss resulting from the impurity/point con- 
tact as the temperature is decreased. This is convenient be- 
cause it allows one to get rid of the contributions proportional 
to L by considering the difference Sjjv ~ Sir, where Suv 
is the entropy in the absence of the impurity/point contact 
(the ultraviolet limit), and Sir is the entropy in the zero- 
temperature limit (the infrared limit). In the IR limit, the point 
contact in the topological theory effectively cuts the edge in 
two. When the edge is cut into two, the topological entan- 
glement entropy Q receives a contribution —2 In T>. Thus 
Sjjv — Sir = In V. As we will see, by definition, T> > 1, so 
this is indeed positive as one expects for a flow caused by a rel- 
evant perturbation. This demonstrates that a "<7-theorem"i2i& 
applies to a point contact across a system in a topological 
phase. Likewise, in quantum impurity problems the impurity 
degrees of freedom end up being screened (at least partially) 
by coupling to the degrees of freedom in the one-dimensional 
bulk. 

These results show that when a (l+l)-d critical system is 
the boundary of a (2+l)-d system in a topological phase, 
the thermodynamic "boundary entropy" of eq. |4] is actually 
the topological entanglement entropy for the subsystems into 
which the system has been dynamically split. This is remi- 
niscent of the generation of entropy by black hole formation. 
(The reader need not worry that by relating a (0+l)-d entropy 
to a (2+l)-d entropy we are forgetting that a boundary cannot 
have a boundary. The point contact is not the boundary of a 
(l+l)-d system, but rather a defect in the middle of it - and, 
therefore, capable of cutting it in two.) 

In fact, the connection between the two kinds of entropy 
goes even deeper. When a topological phase has quasiparti- 
cles with non-Abelian statistics, there is an additional source 
of entanglement entropy. This arises because the Hilbert space 
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for identical particles obeying non-Abelian statistics must be 
multi-dimensional (even when ignoring the position and mo- 
mentum of the particles), so that the system moves around in 
this space as the quasiparticles are braided. For example, in a 
p + ip superconductor, two vortices form a two-state quantum 
system, which therefore have an entanglement entropy ln(2). 
We dub this entropy the bulk entanglement entropy, because it 
comes from the gapped bulk quasiparticles. However, we em- 
phasize that this a subleading piece of the 2+1 -dimensional 
bulk entropy: it is independent of the size of the system, but 
rather depends only on the number and type of quasiparticle 
excitations in a given state. The remarkable connection is that 
this part of 2+1 dimensional bulk entropy can be encoded in 
a 1 + 1 dimensional partition function. In other words, the en- 
tanglement entropy of these bulk quasiparticles is a property 
of the boundary of the system (just like a black hole!). We 
show how to compute such holographic partition functions 
by using conformal field theory. 

In section||l] we discuss the quantum dimensions of a topo- 
logical state and of the various quasiparticle excitations of 
such a state. We give a general way to compute quantum 
dimensions using conformal field theory in section [ni| and 
do this computation for a variety of examples in appendix 
IaI The examples include the Laughlin states for the abelian 
fractional quantum Hall effect, the Moore-Read 17 and Read- 
Rezayi 18 non-abelian quantum Hall states, and the p + ip 
superconductor 19,20 . In section llVl we derive the central results 
of our paper, which define the holographic partition function 
and how it describes the topological entanglement entropy. 
We show in section[V]how these results can be applied to un- 
derstand how a point contact affects the entropy. This com- 
plements recent results of ours, relating point contacts in the 
Moore-Read state and in p + ip superconductors to the Kondo 
problem 2 ^ 2 ^. 



n. QUANTUM DIMENSIONS 

A fundamental characteristic of a topological field theory 
is the quantum dimensions of the excitations. The topolog- 
ical entanglement entropy of the systems studied here on a 
space with smooth boundaries is entirely given in terms of 
these numbers. 

For the quantum Hall effect and related states, one can find 
the quantum dimensions directly from the wavefunctionsSSi. 
In this paper we find it useful to instead compute them us- 
ing the methods of conformal field theory. The formal reason 
why this is possible is that the algebraic structure of rational 
conformal field theory is virtually identical to that of topolog- 
ical field theory (this can be understood precisely by using the 
mathematical language of category theory). A more intuitive 
reason is that the gapless edge modes of the 2+1 -dimensional 
quantum theories studied here form an effectively 1 + 1 dimen- 
sional system. A 1 + 1-dimensional gapless system with linear 
dispersion has conformal invariance, so the powerful meth- 
ods of conformal field theory are applicable. The edge modes 
are in one-to-one correspondence with the modes of the bulk 
system, so the d a computed using conformal field theory are 



those of the bulk system as well. 

We introduce the full conformal field theory formalism in 
the next section. In this section we explain how to define both 
the individual quantum dimensions of the quasiparticles d a , 
and the total quantum dimension I? of a theory. We discuss 
explicitly the simplest non-trivial example, the Ising model. 

The quantum dimension is simple to define. Denote the 
number of linearly-independent states having M quasiparti- 
cles of type a as H a (Af). Then the quantum dimension d a of 
the excitation of type a is given by studying the behavior of 
H a (J\f) for large JV, which behaves as 

H a (JS)ocd% . 

To have non-abelian statistics, one must have d a > 1; these 
degenerate states are the ones which mix with each other un- 
der braiding. In such a situation one typically has d a which 
are not integers. 

As we discussed in the introduction, the topological entan- 
glement entropy is related to the total quantum dimension T>, 
which is defined as 



The sum is over all the types of quasiparticles in the theory. 
Thus to compute T>, one must not only compute the d a but be 
able to classify all the different quasiparticles as well. Since 
da > 1 by definition, for T> to make sense there can only be 
a finite number of different quasiparticles. We must therefore 
confine ourselves to topological field theories where this is 
true. Luckily, the topological theories of most importance in 
condensed-matter physics have this property. 

We will give a simple expression for T> below in terms of 
the modular S matrix of conformal field theory, but here we 
illustrate how it comes about. In the fractional quantum Hall 
effect there is a very natural way of understanding the struc- 
ture of the topological theory. Quasiparticle states are defined 
"modulo an electron". This means that in the topological the- 
ory, two excitations which differ by adding or removing an 
electron are treated as equivalent. Thus the charge-e electron 
is in the identity sector of the theory, since in the topological 
theory the identity and the electron are equivalent. 

A simple example is in the abelian Laughlin states at fill- 
ing fraction v = 1/m. We discuss in depth in the appendix 
how the consequence of modding out by an electron is that 
states differing in charge by any integer times e are identified 
in the topological theory. Only fractionally-charged quasipar- 
ticles result in sectors other than the identity. In the Laughlin 
states, the fundamental quasiparticle has charge —1/m. One 
can obviously consider more than one of these to get larger 
fractional charge, but if one has m of the fundamental quasi- 
particles, they can fuse into a hole (put another way, the hole 
can split apart into rn fundamental quasiparticles). Different 
sectors of the topological theory can therefore be labeled by 
the charges ae/m where a is an integer obeying < a < m. 
There are thus m different quasiparticle sectors. The quantum 
dimension of any quasiparticle in an abelian theory is 1, so we 
have d a = 1 for a = . . . m — 1. Using this with (|5j gives 



4 



the total quantum dimension of a Laughlin state at v = 1/m 
to be 



=l/r 



(6) 



Finding V for non-abelian states requires knowing even 
more about the structure of the theory. Not only does it require 
knowing what the quasiparticles are, but finding the individ- 
ual d a requires knowing what the fusion rules are. The fusion 
rules describe how to treat multi-particle states in the topo- 
logical phase in terms of single-particle states. In the abelian 
case, the fusion rules are simple: for example, two charge- 
e/m particles fuse to effectively give a single charge-2e/m 
particle. The non-abelian structure of the theory arises when 
two particles can fuse in more than one way. 

We illustrate this in the p+ip superconductor. As discussed 
in detail in a number of places 20 - 25 , a vortex has a Majorana 
fermion zero mode in its core. Then two vortices share a Dirac 
fermion (=2 Majorana fermions) zero mode. This zero mode 
can be either empty or filled; thereby, two vortices form a two- 
state system. Denoting the identity sector by /, the vortex by 
a and the filled zero mode by ift, this means that the fusion rule 
for two vortices can be denoted by a ■ a = I+tjj. Two different 
terms on the right-hand-side means that there are two ways to 
fuse two sigma quasiparticles. (More precisely, this means 
that two vortices have a two-dimensional Hilbert space.) One 
finds that the full set of fusion rules here are 



ip ■ ip = I 
a ■ a = I 
ip ■ a = a 



1> 



(7) 



and, of course, fusion with the identity gives the same field 
back. 

The quantum dimensions follow from the fusion rules. Two 
a quasiparticles here can fuse to give either the identity I, or 
the fermion ijj. Thus while H a {l) = 1, these two possibilities 
for fusion of a with itself means that the dimension of the 
Hilbert space for two vortices is H a {2) = 2. Now include a 
third vortex. The dimension of the Hilbert space is not 4, but 
rather is H a (3) = 2. This is because 



a ■ a ■ <j 



(I + ijj) ■ a = a + , 



Fusion is associative, so it makes no difference in which order 
we fuse. Continuing in this fashion, it is easy to see that 



f 2-V/2 

Ha{M) = I 2 (^-l)/2 



AT even 
Af odd 



The quantum dimension of the vortex is therefore d a = \/2. 
Since rp ■ tp = I, H^,(Af) = 1 for all TV, so the quantum 
dimension of the fermion tp is simply d^, — 1. The quantum 
dimension of the identity field is obviously always di = 1, 
so only the vortices exhibit non-abelian statistics. The total 
quantum dimension of the p + ip superconductor is therefore 



V 



p+ip 



l 2 



(V2) 2 



(8) 



As we will discuss in detail in the next section, this all fol- 
lows from studying the Ising conformal field theory, which 
describes the edge modes of this superconductor. 



III. COMPUTING THE QUANTUM DIMENSIONS USING 
CONFORMAL FIELD THEORY 



As discussed in the introduction, the bulk quasiparticles of 
a topological field theory are in one-to-one correspondence 
with the 'primary fields' of a corresponding rational confor- 
mal field theory. Here we explain how to use this correspon- 
dence to give a systematic way to compute the quantum di- 
mensions of the quasiparticles in a topological phase. 



A. Primary fields and quasiparticles 

We start by reviewing some of the conformal field theory 
results we will be using. For an excellent introduction to 
much of this formalism, with detailed applications to the Ising 
model, see Ginsparg's lectures 2 ^. 

We take two-dimensional space to be a disk, so that its one- 
dimensional edge is a circle of radius R. At zero tempera- 
ture the effective two-dimensional spacetime for the confor- 
mal field theory is therefore the surface of a cylinder. We 
describe this cylinder by a periodic variable < 8 < 2tt and a 
(dimensionless) Euclidean time coordinate r. It is often con- 
venient to study the theory on the punctured plane with com- 
plex coordinates z and z, and then use the conformal trans- 
formations z = e T+l6 and ~z — e T ~ 10 to map the results onto 
the cylinder. We will often study the theory at non-zero tem- 
perature 1//?, so that r becomes periodic with period (3. At 
non-zero temperature, spacetime therefore is a torus. All the 
conformal field theories we study are chiral, which means that 
all fields depend only on t + id or z, and not t — iO or z. This 
is possible because of the time-reversal symmetry breaking of 
the 2+1 dimensional theory, coming, for example, from the 
magnetic field required for the Hall effect. 

Many powerful techniques can be used to analyze confor- 
mal field theories. In two spacetime dimensions, conformal 
symmetry has an infinite number of generators. The symme- 
try is generated by the chiral part of the energy-momentum 
tensor T(z) and its antichiral conjugate T(z). The gen- 
erators of chiral conformal transformations are the modes 
of T(z), i.e. the coefficients L n in the Laurent expansion 
T(z) — J2 n L n z~ n ~ 2 on the punctured plane. The energy- 
momentum tensor has dimension 2, so we have normalized 
the modes so that L n has dimension n. The Hamiltonian of 
the non-chiral system on the cylinder is 

H =R {L ° + L ^~6R- 

The constant c is known as the central charge of the confor- 
mal field theory. The term proportional to c arises from the 
conformal transformation of the plane to the cylinder; it can 
be interpreted as the ground-state or Casimir energy of the 
system in finite spatial volume 27rJS2i. Conformal symmetry 
requires that [Lo, L n ] = —nL n , so acting on an energy eigen- 
state with L_„ gives another eigenstate with energy shifted by 
2im/R. Furthermore, an Lq eigenstate has eigenvalue which 
is equal to the dimension of the operator which creates it. It is 
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convenient to define a "chiral Hamiltonian" by 

H = —(L -c/24), 



(9) 



so that H = H + H. 

Since there are an infinite number of symmetry genera- 
tors, the irreducible representations of conformal symmetry 
are infinite-dimensional. Thus one might hope to classify all 
the states of a 1 + 1-dimensional field theory in terms of a finite 
number of irreducible representations of conformal symmetry. 
In other words, there will be a finite number of highest-weight 
states, and all the other states of the theory will be obtained by 
acting on the highest-weight states with the L n with n < 0. 
This is quite similar to decomposing states into irreducible 
representations of any non-abelian symmetry algebra such as, 
for instance, spin. In fact, one often can extend the confor- 
mal symmetry algebra to a larger algebra by also including 
symmetries such as spin, supersymmetry, or even more ex- 
otic possibilities called W algebras. One can then organize 
all the states of the theory into irreducible representations of 
these extended infinite-dimensional symmetry algebras. The- 
ories with a finite number of highest-weight representations 
of such extended symmetry algebras are called rational con- 
formal field theories. The fields creating the highest-weight 
states are called primary fields, and those obtained by acting 
on these with the symmetry generators are called descendants. 

For the theories of interest here, there are in general an infi- 
nite number of primary fields under the conformal symmetry, 
but a finite number under a larger symmetry group. There is 
a marvelous physical interpretation of the presence of this ex- 
tended symmetry algebra. We discussed in the previous sec- 
tion how in the Hall effect, the topological field theory arises 
by considering all the states "modulo an electron". In the 
edge conformal field theory, the symmetry algebra can be ex- 
tended by including the electron annihilation/creation opera- 
tors. Thus the primary fields of the conformal field theory are 
in one-to-one correspondence with the quasiparticles in the 
topological phase. Descendant states arise by attaching elec- 
trons or holes to the quasiparticles - acting with the electron 
creating or annihilation operator moves one around inside an 
irreducible representation of the extended symmetry algebra. 
We will give explicit examples below of how this works. 



B. Quantum dimensions from the fusion rules 

There are several equivalent ways of using conformal field 
theory to compute the quantum dimensions of the quasiparti- 
cles. All amount to finding the fusion coefficients. Since all 
the fields of the theory are expressed by acting with the sym- 
metry generators on the primary fields, the operator-product 
expansion of two primary fields can be written as a sum over 
the primary fields. Namely, the "fusion rules" of a conformal 
field theory are written as 



E 



(10) 



where the fusion coefficients N° b are non-negative integers, 
which count the number of times the primary field cf> c appears 
in the operator-product expansion of <fr a and <fi b . This algebra 
is associative 



and one has N£ b = 



ba 



N^. For all the 



theories considered here, we have N£ b = or 1, but it is 
possible to have larger integers in more complicated theories. 

Each quasiparticle in the topological field theory corre- 
sponds to a primary field, and the fusion of these quasiparti- 
cles is the same as that of the corresponding primary fields. 
Thus to have non-abelian statistics, N£ b for some a and b 
must be non-zero for more than one c. A simple example 
of non-abelian statistics is given by the Ising model, which 
has three primary fields: the identity field /, the fermion 
ip, and the spin field a which have the fusion rules given 
in (|8j above. The Ising CFT describes the edge excitations 
of a p + ip superconductor, with the vortices correspond- 
ing to (7i22i££ This can be shown, for example, by using the 
Bogoliubov-de Gennes equations for the superconductor, as 
reviewed in ref. The non-zero fusion coefficients are 
N* a = N* c = N^, = = 1. These fusion rules are 
compatible with the spin-flip symmetry a — > — a, ip — > ip. 

We detailed in the previous section how to compute the 
quantum dimensions for the Ising fusion rules. The general 
procedure is similar. If we fuse M <p a fields together and use 
dint recursively on the right-hand side: 

0a • 0a " • • ■ " 0a = K'aKl, ' ' ' K^K^ ( U ) 



This is the product of M — 1 copies of the matrix 

(Qa)t = Kb ■ 



(12) 



In the M — ► oo limit, the product will be dominated by the 
largest eigenvalue of Q a . This is the quantum dimension d a . 

This eigenvalue can easily be written in terms of the mod- 
ular S matrix, which we will define and discuss in the next 
subsection lHICI A profound result of conformal field theory 
is the Verlinde formula^^, which expresses the fusion rules 
in terms of the elements of S, namely, 



ab 



S 3 a S J b Sj 







<? 3 
o 



(13) 



where denotes the identity field and the sum is over all pri- 
maries. In general S is unitary and hermitian; in all of our 
examples S is also real and hence symmetric. We therefore 
give the formulas for this case; they are simple to generalize. 

To find the eigenvalues and eigenvectors of Q a , we multiply 
the Verlinde formula by S% / Sq and sum over b. This yields 



.S 



qk cc 

0„ lJl 



b 



i r\ \ c a k 

\yta)b gk ~ Qk qc 



I] 



ck Oi 
o O 



The 6th element of the eigenvector of Q a with eigenvalue 
S% I Sq is therefore given by S^ /Sq. The modular S matrix is 
the matrix of the eigenvectors of Q a . The largest eigenvalue 
is the one with k = 0, so we have 



CO 

d = -2. 



(14) 
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There is a simple expression for the total quantum dimen- 
sion T> in terms of the modular S matrix. Using (I14> along 
with the fact that S is unitary and symmetric gives 




This formula will prove useful in relating the topological and 
thermodynamic entropies. 

C. The modular S matrix 

Rational conformal field theories have a variety of profound 
mathematical properties. We have exploited one of them in 
using the Verlinde formula Jl 31 to give a compact expression 
(114-i for the quantum dimensions in terms of the modular S 
matrix. In this section we will discuss what S is and how to 
compute it. 

To define S and to use the Verlinde formula, we study con- 
formal field theory on a torus, i.e. in finite size in both the 
Euclidean time and spatial directions. Physically, this corre- 
sponds to a non-zero temperature, so that the Euclidean time 
coordinate is periodic with period (3. One must compute the 
partition function in the sector of states associated with each 
primary field a and its descendants. In mathematical lan- 
guage, this is called a character, and is defined as 

Xa(q) = tx a e~ m = q- c/24 \x a q L \ (16) 

where q = e - 27 W^ ji j s the spatial length of the system, and 
the trace tr a is over all the states in the irreducible represen- 
tation of the extended symmetry algebra corresponding to the 
primary field a (i.e. the highest-weight state and its descen- 
dants). The partition function of the chiral theory is then of 
the form 

Z = Y J N a Xa(q), (17) 

a 

where the N a are integers, representing how many copies of 
each primary field appears in a given theory. In a non-chiral 
theory, one defines analogous antichiral partition functions 
X a (q), and the full partition function of a rational confor- 
mal field theory is given by a sum over products of the form 
Xa (q)Xb (?)> a g a in with integer coefficients. In both chiral and 
non-chiral rational conformal field theories, the sums are over 
a finite number of characters. We will give explicit examples 
of these characters below. 

The modular S matrix describes how the characters behave 
when one exchanges the roles of space and time: 

Xa(q) = Y,S b aX b(q). (18) 

b 

where q = e~ 27TR /P. Keep in mind, however, that a 1+1- 
dimensional quantum theory in Euclidean time is equivalent 
to a two-dimensional classical theory. Going back from the 



two-dimensional classical theory to a 1+1-dimensional quan- 
tum theory, one is free to choose either of the directions to be 
space and the other to be time. We thus conclude that the the- 
ory should be invariant under modular transformations of the 
torus. Hence, 

Z = Y J NaXa(q)- (19) 

a 

However, the N a = ^ b S^Ni, are usually not integers; we 
will have much more to say about this in section llvl 

Computing the characters is a problem in the representation 
theory of the extended symmetry algebra. In the cases of in- 
terest here the computation is straightforward; we give several 
examples in appendix|X| Finding the integers N a for a given 
bulk topological state is the crucial computation in this paper, 
and we discuss how these arise, and what they have to do with 
the entropy, in section lTvl 

IV. HOLOGRAPHIC PARTITION FUNCTIONS 

In this section we arrive at a central result of this paper. We 
first define the edge entropy in terms of the chiral conformal 
field theory describing the edge modes. The edge entropy is 
not merely reminiscent of the topological entanglement en- 
tropy Q discussed in the introduction, but has an identical 
universal part — ln(2?). We then extend the correspondence 
between entanglement entropy and thermodynamic entropy 
by studying the bulk entanglement entropy, which arises from 
the different fusion channels of the bulk quasiparticles. We 
show that both of these entropies can be encoded in a sin- 
gle holographic partition function, which we interpret as the 
topological entanglement entropy of a region of a system in a 
topological phase. 



A. Edge and "boundary" entropies 

We are studying 2+1 dimensional systems which are 
gapped in the bulk but gapless at the edge. Their edge modes 
are described by chiral rational conformal field theories. In 
the previous sections, we have defined and discussed the char- 
acters of conformal field theories, which are essentially chiral 
partition functions. Here we use these to characters to define 
the edge entropy. 

To motivate our definition of the entropy, we first discuss 
a closely-related system, a non-chiral rational conformal field 
theory with boundaries. This means that spacetime at zero 
temperature is a strip of width ttR instead of a cylinder. 
At non-zero temperature, spacetime becomes a finite-width 
cylinder with period f3 in the Euclidean time direction, instead 
of a torus. The left and right movers on the strip are coupled 
by the boundary conditions, as illustrated schematically in fig- 
ure^ This means the chiral and antichiral conformal symme- 
tries are also coupled. Conformal invariance is preserved only 
with certain boundary conditions, and even then, only a sin- 
gle algebra remains^. This symmetry algebra is exactly the 
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FIG. 1 : The edge modes with spacetime a cylinder (illustrated on the 
left) have only one chirality. The physics of these is closely related 
to that of a non-chiral system with spacetime a strip (illustrated on 
the right), where the boundary conditions couple the left and right 
movers. 



same as that of a single chiral theory. The partition function 
for the model on the strip is then a sum over characters of 
a single algebra. These characters are exactly those in dl 6I >. 
and the partition function takes the same form dl7> . with N a 

» 12 29 

tntegesiaifii. 

The connections between the chiral theory on the cylinder 
and the non-chiral theory on the strip go even deeper. The 
boundary conditions which preserve conformal invariance are 
in one-to-one correspondence with the primary fields of this 
chiral rational conformal field theory, so each can be labeled 
by the same indices we use to label primary fields. (More pre- 
cisely, conformal boundary conditions form a vector space, 
basis vectors of which can be labeled by the primary fields.) 
Then a key result of Cardy's is that with appropriate regular- 
ization, the partition function Zjk on the strip with boundary 
conditions j and k at the two ends isi^ 

Z 3 k=Y. N 0kXa{q) (20) 
a 

where N? k is the same Nj k which appears in the fusion rules! 
Furthermore, by using the Verlinde formula it is simple to find 
boundary conditions which result in a partition function given 
by a single character \a- 

These results were used by Affleck and Ludwig to define 
and compute what they called the "boundary entropy'iii Con- 
sider an RCFT on the finite-width cylinder with boundary con- 
ditions j and k, so that the partition function is Zjk. Then take 
the width of the cylinder R — ► oo while holding the temper- 
ature, and hence the radius (3 of the cylinder, fixed at some 
non-zero finite value. Then, on general grounds, one must 
have 

lnZ jk = fR + ln(g) + ... (21) 

where we neglect terms that vanish as R — > oo. This is pre- 
cisely the form (|4}. We can compute / and g by using a mod- 
ular transformation, because the limit R/f3 — > oo corresponds 



to q = 0. From Jl 8i we have 

a , b 

By definition, in this limit, the character Xa behaves as 

So unless j and k are such that ^ Q Nj k S® = 0, the identity 
character x dominates and one has 

a 

It is then natural to interpret g as a ground-state degeneracy, 
and ln(g) as a boundary entropy, since clearly this subleading 
term in the full entropy depends on the boundary conditions 
j and k. However, g is not necessarily an integer, so as we 
discussed in the introduction, it cannot, strictly speaking, be a 
ground-state degeneracy. Neither is it solely associated with 
the boundaries of the strip, since there is no particular reason 
that subleading terms from the 1 + 1 dimensional bulk cannot 
contribute to ln(<7). 

We now return to studying topological theories whose edge 
modes are described by chiral conformal field theories on the 
torus. Although the torus has no boundaries, there are still 
boundary conditions on the fields. For example, in the p + ip 
superconductor, the edge fermions have antiperiodic or peri- 
odic boundary conditions around the spatial cycle of the torus 
depending on whether there are an even or an odd number of 
vortices in the two-dimensional bulk. So let us first consider 
the case where there are no bulk quasiparticles. In the lan- 
guage of topological field theory, this corresponds to trivial 
topological charge on the disk. The partition function of the 
edge RCFT is then 

z o = Xo(o) 

where, as above, the label on \o means the identity sector. 
To extract the piece of interest from this, we take the same 
R//3 — > oo limit. Z Q then can be expanded in the form d2"TT i. 
i.e. 

m(Zo) = /ofl+ln(7o), (23) 

where 7 in this chiral theory is the analog of the g in J2H in the 
impurity theory. Following through the same modular trans- 
formation described in the preceding paragraph, we obtain the 
subleading contribution to ln(x ) in this limit: 

7o = AM = S° . 

Recall, however, that we showed in {T3) that £g = l/V, 
where V is the quantum dimension we have gone to great 
lengths to compute. We thus conclude 

S edge ^ln(7 ) = -ln(2?) (24) 
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as the universal part of the edge entropy. There is no particular 
reason for 70 to be an integer, since it is a subleading term in 
the large R//3 expansion. 

Note that 5 e( jg e is identical to the universal piece of the 
topological entanglement entropy in 0. Although the two 
quantities are defined in different ways, this is hardly a coin- 
cidence. 5 e( jg e is defined for a system with a real edge, while 
the topological entanglement entropy is defined on the bound- 
ary between two regions A and B, not a physical edge. Nev- 
ertheless, one expects the two to be the same. Edge modes 
arise formally in a topological theory to cancel the chiral 
anomaly 30 . The topological entanglement entropy arises for- 
mally by integrating out the system beyond the boundary, i.e. 
subsystem B. When integrating out the degrees of freedom 
of B, one of course integrates out an anomaly-free theory, so 
this must include the edge modes of B as well. The remaining 
A theory then, in some sense, must include "edge modes" as 
well which are required to cancel those on B. This required 
cancellation is one way of seeing that systems A and B are 
indeed entangled. Thus, the value of the topological entan- 
glement entropy should be identical to our computation of the 
edge entropy. We have shown by direct computation that this 
is indeed the case. 



B. Bulk entanglement entropy 

The correspondence between thermodynamic and entangle- 
ment entropies goes much deeper. Let us now consider the 
bulk entanglement entropy. This arises in non-abelian topo- 
logical phases because of multiple possible fusion channels. 
For example, two vortices in a p + ip superconductor form 
a two-state quantum system, so their entanglement entropy is 
ln(2). 

We can derive a simple formula for the bulk entanglement 
entropy in general. We define the topological degeneracy of a 
given state to be the number of ways one can fuse the quasipar- 
ticles present to get a particular overall fusion channel. This is 
easily written in terms of the matrix Q a defined in d!2i . When 
there are n a of each of the quasiparticles labeled by a, define 
the matrix T by 

T = l[(Qa) n " (25) 

a 

It makes no difference in which order we multiply the Q a be- 
cause different Q a commute with each other (as we showed 
from the Verlinde formula, they have the same eigenvectors). 
The topological degeneracy of channel c is then the c, entry 
of this matrix T. Note that this is a simple generalization of 
the way we computed the quantum dimension to the case in 
which different types of particles are allowed to be present. 

Let us illustrate this in terms of the p + ip superconduc- 
tor, where there are three states in the topological field the- 
ory, labeled by /, ip and a. If there are no bulk quasiparticles 
present, we have degeneracy 1 in the identity channel. If we 
have only ip quasiparticles, the fusing is simple: for an even 
number of ip quasiparticles the degeneracy is 1 in the iden- 
tity channel, while for an odd number the degeneracy is 1 in 



the ip channel. It gets more interesting for bulk a quasiparti- 
cles. A single bulk a quasiparticle corresponds to degeneracy 
1 in the a channel. However, two bulk a quasiparticles does 
not correspond to degeneracy 2. Since a ■ a = I + ip, two 
bulk a quasiparticles corresponds to degeneracy 1 in both the 
/ channel and in the ip channel. Degeneracy 2 in the a chan- 
nel occurs for three bulk a quasiparticles. In general, for an 
even number 2M of bulk a quasiparticles, we have degener- 
acy 2 M_1 in the I and tp channels while for an odd number 
2M + 1, we have degeneracy 2 M in the a channel. Including 
any number of ip quasiparticles doesn't change the degenera- 
cies (except for M = 0). 

The bulk entanglement entropy arises from the uncertainty 
in knowing which quantum state the bulk quasiparticles are 
in. For example, if there are two bulk a quasiparticles, we do 
not know without doing a measurement whether they are in 
the / or ip channel. The entropy is therefore ln(2). But what 
if there is just one a quasiparticle in a region? It can be en- 
tangled with a a quasiparticle in another region of the droplet, 
giving a total entropy of ln(2). Since entropy is additive, the 
only consistent possibility is for the single a quasiparticle to 
have entropy ln(\/2). In general, a single quasiparticle of type 
a has entropy ln(d ), where d a is the quantum dimension dis- 
cussed at great length above. The general expression for the 
bulk entropy can be expressed easily in terms of the topologi- 
cal degeneracies, namely 

5 bulk = ln (E T °c4). (26) 

This can be simplified by rewriting d c and T c q in terms of the 
modular S matrix. We can use the Verlinde formula Jl 3i for 
Q and the relation d c = S^/Sq. The latter is the eigenvector 
of any Q a , having eigenvalue S®/Sq = d a . Thus it is an 
eigenvector of T as well, so 

5 bulk = E"aln(d a ) • (27) 

a 

Each quasiparticle of quantum dimension d a contributes 
ln(d ) to the bulk entanglement entropy. Note that this "bulk" 
entropy is extensive not in the size of the sample, but rather in 
the number of bulk quasiparticles. 

C. The bulk and edge entropies from the holographic partition 
functions 

We show here that the entire topological entanglement en- 
tropy can be seen in terms of the edge conformal field theory. 
We interpret this as an example of holography. In a topolog- 
ical field theory, the degrees of freedom live on the edge. In 
the fractional quantum Hall effect, these are, of course, the 
edge modes we have been discussing all along. Holography 
simply means that all the information of the bulk topological 
field theory is encoded in the edge modes. 

As we have discussed above, adding or removing bulk 
quasiparticles changes the boundary conditions on the edge 
modes. This will change the chiral partition function for the 
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edge. To illustrate this, consider the p + ip superconductor, 
where there are three characters: xi< Xv anc ' Xcrt w i m a mod- 
ular S matrix given in ED. The presence of a quasiparticles 
(the vortices) in the bulk has an important effect on the edge. 
As we reviewed in ref. |22j, an odd number of a quasiparti- 
cles changes the boundary conditions of the edge Majorana 
fermion from antiperiodic around the circle to periodic. For 
a fermion with periodic boundary conditions on the 1 + 1 di- 
mensional torus, the chiral partition function is \ a . Thus we 
define the holographic partition function for a single bulk u 
quasiparticle to be Z a = Xa- This partition function yields 
the correct entanglement entropy. Taking the L/fi^oo limit 
as in d23l > yields 



For ap + ip superconductor with one vortex in region A, the 
entanglement entropy is indeed 

5 bulk + 5 edge = MV2) - ln(2) = - \n(V2). 

We define the holographic partition function for a single bulk 
ip quasiparticle in a similar fashion: = Xi>- The corre- 
sponding entropy is then ln(Sf) — ln(l/2) = — In 27. This 
is just <S ec |g e , which is indeed what we expect, because there 
is no bulk topological entropy for ip quasiparticles. 

Given our definition of the bulk entanglement entropy, it 
should be obvious how to define the holographic partition 
function for an arbitrary number of bulk particles. The ef- 
fect of adding a bulk quasiparticle is given by the same fusion 
rules as for the primary fields. For example, in the p + ip su- 
perconductor, the fusion rule ip ■ ip = I means that with no a 
quasiparticles and an even (odd) number of ip quasiparticles, 
we define the holographic partition function to be xi (Xi/>)- ^ n 
both cases the entropy remains — In V, since there is no bulk 
entanglement entropy. A less trivial case is when we have two 
bulk a quasiparticles. Here the holographic partition function 
follows from the fusion rule a ■ a = I + ip, so we get 

Zla = Xl + Xi>- 

For three, the topological degeneracy is 2 in the a channel, so 
we get 

Z-ia = 2x CT , 

and so on. 

The definition of the holographic partition function for ar- 
bitrary numbers of quasiparticles in the bulk should now be 
obvious. In general, for n a bulk quasiparticles of each type a, 
the holographic partition function is 

Z = Y, T 0aXa(q) (28) 

a 

where T is the topological degeneracy matrix defined in d25l >. 
This is the partition function obtained by computing the path 
integral for the appropriate 2+1 dimensional topological field 
theory (e.g. Chern-Simons theory) with spacetime D x Si, 



where D is a disk of radius R, and Si a circle of radius j3. To 
get the correct degeneracies T$ a one needs to insert n a Wilson 
loops of each type a which puncture the disk and wrap around 
the Si. 

As advertised, the holographic partition function d28i gives 
the correct entanglement entropy in the the R//3 — > oo limit. 
In this limit, 

Z = J2ToaS Q a X (q) + ... . 

a 

We define the total thermodynamic entropy S via our usual 
expansion 

In(Z) = fR + S + ... . 

Using db = Sq/Sq along with the expressions for the bulk 
and edge entropies j26t and &27\ gives 

S = ln(^T 0a S o °) 

a 

= ln(^T 0a rf a )+ln(S °) 

a 

= ^n a ln(d Q )-m(2?) . (29) 

a 

Thus we see that indeed once we have taken the R//3 — > oo 
limit, we have 

S = S bu\k + 5 edge ■ ( 3 °) 

The holographic partition includes both bulk and edge entan- 
glement entropies. 

We close this section by returning to where it started, a dis- 
cussion of the relation between chiral conformal field theory 
with spacetime a torus and non-chiral conformal field theory 
on a finite-width cylinder. The partition functions in both 
cases are expressed in the form 

Z = Y,NaX« 

a 

where the N a are integers, describing how many times each 
character appears. In the latter case, however, there is one 
important restriction which we need not demand of the holo- 
graphic partition function. For a partition function describing 
a unitary field theory on the finite-width strip with physical 
boundary conditions, the identity character usually appears at 
most once (i.e. Ao = or 1). The reason is that there is 
only one identity field in a unitary theory, so only with pe- 
culiar boundary conditions can one have Nq > 1. In the 
holographic partition function, Nq can be any non-negative 
integer: it is the number of ways the bulk quasiparticles in the 
2+1 dimensional topological theory can fuse to get the iden- 
tity channel. Thus although the partition functions are defined 
in a very similar fashion, they really are different objects. 

V. DYNAMICAL ENTROPY LOSS AT A POINT CONTACT 

The relation between the thermodynamic and entanglement 
entropies manifests itself most dramatically through the dy- 
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namics of a point contact connecting two points on the bound- 
ary of the disk. Quasiparticles can tunnel between the two 
points connected by the contact, thereby perturbing the edge 
modes at these two points. This can change the edge en- 
tropy. When the tunneling operator is a relevant perturbation, 
it causes the system to flow to an infrared fixed point at which 
any edge excitation which is incident upon the contact neces- 
sarily tunnels across the system. In other words, an imaginary 
boundary between the two halves evolves dynamically into a 
real edge. As we shall see below, this process causes the en- 
tanglement entropy between the two halves to be, in a sense, 
"reborn" as the actual thermodynamic entropy of the two re- 
sulting droplets. We have discussed a point contact in a non- 
abelian topological state in depth in a companion paper—, and 
the results in this section complement those results. 

As mentioned in the introduction, if we compute the den- 
sity matrix for half of a system by tracing out the degrees of 
freedom of the other half, the resulting density matrix has an 
entanglement entropy 

Sa = o,L — In T> + . . . 

The second, universal term is negative. What it tells us is that 
there is a little less uncertainty in the state of half of the sys- 
tem than one might have naively expected purely from local 
correlations 3 ' 4 . A priori, this has absolutely nothing to do with 
the thermodynamic entropy of the system since it is obtained 
from an arbitrary division of the ground state wavefunction 
alone, with nothing about the spectrum taken into account. 
However, when the strong tunneling fixed point is reached, 
the disk is effectively split into two. The imaginary division 
into two halves becomes a real division. Entropy is lost as a 
result of this process because the edge entropy for two disks 
is —2 ln(Z>), compared to — In (2?) for a single disk. As in the 
case of the entanglement entropy, this extra negative contribu- 
tion is due to the fact that we know more about the system: 
the negative contribution — ln(2?) to the entropy of a disk rep- 
resents the information that we know about the system when 
we know that it has trivial total topological charge. However, 
one half of the system could have arbitrary topological charge 
a, so long as the other half has the compensating topologi- 
cal charge a. If the point contact breaks the system into two 
systems, each with trivial topological charge, then this uncer- 
tainty is lost. The information which we now know about the 
two disks is —2 In(X>). 

Defining Sjjv to be the entropy in the limit in which there 
is no tunneling across the sample (the UV limit), and Sir to 
be the entropy in the limit in which tunneling is strong (the IR 
limit), we have 

Suv - Sir = - ln(2?) - (-2 ln(D)) = m(X>) . (31) 

Here, we have assumed that there are no quasiparticles in the 
bulk, and that the system breaks into two halves each with no 
quasiparticles. However, the entropy change is the same even 
when quasiparticles are present, which we will discuss below. 

Let us first illustrate entropy loss due to a point contact in 
the abelian case, when the edge theory consists of a free bo- 
son. This has been widely studied in the context of the Laugh- 



lin fractional quantum Hall states . By using a series of map- 
pings, the chiral problem on the torus can be shown to be ex- 
actly equivalent to a non-chiral problem on the finite-width 
cylinder. In other words, the analogy we made in the previous 
section is an exact equivalence in this case. The point contact 
results in a boundary with Neumann boundary conditions on 
the boson in the UV limit, and Dirichlet boundary conditions 
in the IR 32 . The edge entropies can be computed in both cases, 
and one obtains 32 - 33 for v — 1 /m 

S uv - Sir = ln(Vm). 

This is exactly what one obtains from our results by assuming 
that the IR limit consists of two separate disks. We show in 
iA6\ that T> = *Jm for the v = 1/m Laughlin states. Thus 
we indeed have Suv — Sir = ln(Z>) as in d3 II . 

Using the results derived above, we can check d3 1 1 in sev- 
eral non-abelian cases. As we showed in ref.|23, a point con- 
tact in the p + ip superconductor is equivalent to the single- 
channel anisotropic Kondo problem. The entropy change is 
well-known here. The Kondo problem consists of a single im- 
purity spin antiferromagnetically coupled to an electron gas. 
In the single-channel spin- 1/2 case, the spin is decoupled in 
the UV limit, while it is completely screened in the IR limit. 
Therefore Suv — m (2), while Sir — ln(l) = 0. We showed 
in (|8j that for the Ising model, 2? P + Jp = 2, so i31\ holds. 

A point contact in the Moore-Read Pfaffian state at v = 1/2 
(or v — 5/2) turns out to be a variant of the two-channel 
overscreened Kondo problem. In ref. |2l|, we showed that the 
entropy loss here is Suv — Sir = 21n(v2). In ( IA10t . we 
show that T>mr — 2\/2, so again (13 1 1 holds. For a point 
contact in a Read-Rezayi state or for SU(2)k for general k, 
it is not yet known how to map the tunneling operator onto a 
Kondo-like problem, so we do not have an independent result 
against which to check the prediction which can be obtained 
from i31\ along with the quantum dimension in JA13i . 

The entropy change J3 1 1 can also be seen in terms of the 
holographic partition functions. The partition function for two 
independent systems is simply the product of the partition 
functions for the two systems. Thus when the point contact 
splits the system into two, the holographic partition function 
should be the product of two holographic partition functions. 
Let us first examine this in the case where there are no bulk 
quasiparticles. In the IR limit, the system decouples into two 
independent discs, each with no bulk quasiparticles. Thus the 
partition function is Zq = \o (<?) for the UV, and 

z o,o = XoMXoM 

for the IR. We have denoted qi = exp(— 2irj3/Ri) and q-i = 
cxp(— 2irf3/B,2), where R± and R2 are the sizes of the discs 
into which the point contact splits the system in the IR limit. 
The universal piece of the entropy extracted from Xa{q) is 
independent of how q is taken to 1, so we obtain entropies of 
— ln(2?) for the entropy associated with each of the two disks. 
Thus from Zo,o we obtain Sir = —2 ln(Z>), as required. 

When there are bulk quasiparticles, they must end up on 
one side or the other once the disk is split. The UV holo- 
graphic partition function for a quasiparticle of type a is de- 
noted Xa(q), which from d29t gives the entanglement entropy 
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FIG. 2: In the UV, bulk quasiparticles can be anywhere. In the IR, 
they must end up on one of the two sides. For two bulk quasiparticles 
of type a, they can both end up on one side, or on separate sides. 

of Sjjv = ln(^a) — ln(X>). When there is one bulk quasiparti- 
cle, it must end up on one half or another, but it cannot end in 
both. There are therefore two possibilities for the IR partition 
function: 

Z a ,0 = XaOOXofe)) or Z 0,a = Xo(9i)Xafe). 

We cannot predict which one is realized, but one or the other 
will happen. Either of the two yields the correct entanglement 
entropy 

S IR = la(d a )-2hx(V). 

For multiple particles, we must apply the fusion rules just 
as we did for the UV holographic partition functions. For two 
quasiparticles of type a, there are three possibilities for how 
they might end up, as illustrated schematically in figure M 
If one quasiparticle ends up on each disk, then the IR holo- 
graphic partition is Z a . a = Xa(qi)Xa(q2)- If both quasiparti- 
cles end up on the same disk, then e.g. 

Z 2 afi = Z 2a (qi)Z (q 2 ) = ^2 N aaXb{qi)^j XofeO- 

For example, for the p + ip superconductor with two bulk a 
quasiparticles, we have 

Z2a,o = {Xl(qi) +Xv(9i))X-r(<72), 

Z a ,„ = X<x(<?l)X<x(<?2)- 

Thus the general rule is simple: for a given subdivision, the 
holographic partition on each half is the same as it would be 
for an isolated disc. Then multiply the characters for the two 
separated halves together. 

With this holographic partition for the IR, the bulk part of 
the entropy remains ^ Q n a ln(d a ) in the IR. It does not mat- 
ter which half the bulk quasiparticle is on: it still contributes 



ln(d a ) in accord with the intuition that a point contact should 
not affect the bulk of the sample. As we have seen in ( 13 II . the 
edge entropy changes, doubling from — ln(Z>) to —2 ln(Z>). 



VI. DISCUSSION 

The ideas presented here grew out of our attempts to 
understand quasiparticle tunneling at a point contact in a 
Moore-Read Pfaffian non-Abelian quantum Hall state^S. A 
key technical challenge in this problem was to formulate a 
bosonized description of the tunneling of twist fields (or Ising 
spin fields) a between Majorana fermion edge modes. We 
found that this could be done if one introduced a spin- 1/2 
degree of freedom at the point contact which served a 'book- 
keeping' purpose of keeping track of how many a fields had 
tunneled. We thereby saw that the tunneling Hamiltonian for 
this problem could be mapped onto various versions (depend- 
ing on filling fraction of the quantum Hall state) of the Kondo 
model. Therefore, the issue of entropy loss due to the point 
contact was inescapable, although the same issue arises in 
simpler cases such as Abelian quantum Hall states but could 
be ignored there. Our interpretation 21 , on which we have elab- 
orated here, is that this entropy loss cannot be ascribed to 
the point contact per se. Rather, it should be understood as 
a (2+1) -dimensional entanglement entropy. When tunneling 
at a point contact causes a Hall droplet to be effectively bro- 
ken in two, the topological entanglement entropy between the 
two sides becomes actual thermodynamic entropy. Since this 
entropy is negative, S = — In T>, the entropy of the system 
decreases. 

One implication of our results is that an attempt to for- 
mulate a tractable representation of tunneling in other non- 
Abelian quantum Hall states, such as the Read-Rezayi states, 
will require introducing an 'impurity' degree of freedom at the 
point contact, analogous to the Kondo spin in the Moore-Read 
Pfaffian case, to give us the correct entropy loss. 

It would be interesting to see which of these results apply 
to gapless 2+1 dimensional theories in a topological phase. 
Such theories are often closely related to conformal field the- 
ories, but these theories are not describing the edge modes, 
but rather the behavior in two-dimensional space. The en- 
tanglement entropy for such systems 34 is quite similar to the 
entropies discussed here, and it would be interesting to pursue 
the analogies further. 

It would also be interesting to explore the connection of 
these results with those motivated by gauge theory, string the- 
ory, and black holes. Results quite analogous to ours have 
been derived using the AdS/CFT correspondence 35 . The con- 
formal field theories here are simpler than those which gener- 
ally arise from the AdS/CFT correspondence (ours are unitary 
and rational), and it would be interesting to understand how 
quantum Hall physics fits into this more general setting. 
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APPENDIX A: EXAMPLES OF QUANTUM DIMENSIONS 
1. Ising model/p + ip superconductor 



for integer n. The operator e l7V appears frequently in the 
Coulomb-gas approach to 2d statistical mechanical models, 
and in string theory as a vertex operator. 

With the conventional normalization of tp, Vy has scaling 
dimension j 2 /2. This means that the eigenvalue of H acting 
on the state created by V-y is 2n('j 2 /2 - 1 /24) / R. More stan- 
dard calculations show that the operator product expansion is 




The Ising conformal field theory is equivalent to a free 
massless Majorana fermion. The partition function on the 
torus - and hence the characters of interest - can be computed 
using the properties of free fermions; the main subtlety is that 
to compute characters for all the primary fields I,ip,a one 
needs to utilize both periodic and antiperiodic boundary con- 
ditions on the fermion. The connection between the different 
characters and the different boundary conditions is reviewed 
in depth in, for instance, ref. |3 In the P + ip superconduc- 
tor, inserting a a quasiparticle (i.e. a vortex) in the bulk flips 
the boundary conditions on the fermion between periodic and 
antiperiodic around the disk (see e.g. ref. 22). 

It is natural and convenient to write conformal-field-theory 
characters in terms of Jacobi theta functions. The behavior of 
theta functions under modular transformations is simple, so 
it is then straightforward to work out the modular S matrix. 
Explicit expressions for the characters can be found in Refs. 
26 36. The modular S matrix for the Ising model is2& 



(Al) 



where we write the fields in the order /, ip, a. Note that S 2 = 
1, as required. 

Using H41 . we find dj = 1, d^, = 1 and d a = y2, as we 
saw directly from the fusion rules. Using J15I . we find the 
total quantum dimension T) p+ i V = 2, as we saw in 



2. Free boson/Laughlin states 

Before studying the non-abelian Moore-Read Pfaffian state, 
it is first useful to study the abelian Laughlin state at v = 1 /3. 
The latter is the first state in the a sequence of quantum Hall 
states, with the Moore-Read Pfaffian second. The quasipar- 
ticles and quasiholes in the Laughlin state at v = 1/3 are 
abelian anyons, with charges ±1/3. 

The conformal field theory describing the edge modes of 
the Laughlin states is simply a free chiral boson <p(z)2m. Here 
(and in most situations of interest), the boson is compact, 
meaning in our chiral context that operators must be invari- 
ant under the shifting ip(z) to ip(z) + 2nr for some "radius" 
r (not to be confused with the radius R of the disk; calling 
r a radius is a relic of string theory). Then standard manipu- 
lations show that a free boson has central charge c = 1, and 
that the primary fields on the punctured plane under conformal 



symmetry are 



26 



Vy(z) 



j7V>(z) 



n 
r 



< ■■■■-■( ' ■ -i = (z — w) '"I 



7<5 i~{ip(z)-\-i8<p(w) 



(A3) 



To obtain correlation functions on the cylinder, it is usually 
most convenient to first work them out on the punctured plane, 
and then conformally transform them to the cylinder. In gen- 
eral, a primary field ~z) of dimensions (h, h) transforms 
under the conformal transformation z' — f(z) as (again, see 



e.g^) 



df 

dz 



h 



df 

dz 



$(f(z),f(z)) 



For the transformation from the punctured plane to the cylin- 
der, we have r + ia = R\h(z)/2tt, so for the operators V 1 {z), 
we have 



R 



7 2 /2 



(A4) 



(A2) 



The fact that there is only one primary field on the right- 
hand-side of the operator product expansion jA3> means that 
the quantum dimension of all the operators are d 7 = 1. 
Computing the total quantum dimension T> for a free boson is 
more work, because there are an infinite number of primaries 
if we use only conformal symmetry since the operator product 
expansion for N operators V 1 gives the operator V/v 7 . One in- 
teresting symmetry generator is the U(l) current d z ip. In the 
fractional quantum Hall effect, this is the electrical charge cur- 
rent. To define T>, we must find the largest possible symmetry 
algebra. For any rational r 2 , one can extend the symmetry 
algebra by some integer-dimension operator and obtain a fi- 
nite number of primary fields, thus giving a rational conformal 
field theory. 

However, for the v — 1/ra Laughlin states in the frac- 
tional quantum Hall effect, the extended symmetry is even 
larger. As we discussed above, one needs to extend the al- 
gebra by the electron annihilation/creation operators. These 
have dimension m/22£. A field with half-integer dimension is 
fermionic; symmetries with such generators occur for exam- 
ple in the Gross-Neveu model for an odd number of Majorana 
fermions 37 . (Incidentally, in these Gross-Neveu models, the 
massive quasiparticles end up having a quantum dimension of 
V53) For the simplest non-trivial case m — 3, the generators 
are of dimension 3/2. The symmetry here is supersymmetry^. 
It is called Af — 2 supersymmetry, because there are two gen- 
erators: the electron creation and annihilation operators. The 
entire sequence of Read-Rezayi states 18 , of which the v = 1/3 
Laughlin state is the first member and the Moore-Read state 
the second, has Af = 2 supersymmetry as its extended sym- 
metry. 
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Using the operator product expansion, it is simple to see 
what the primary fields are for the v = 1/3 Laughlin state. 
The supersymmetry generators are G ± = e ±l ^ v , which in- 
deed have dimension 3/2 and charge ±1. From ( IA3I . we have 



G ± (z)V r (w) = (z-w) ±V *W 1±V s(w)(l + ...) (A5) 

where the terms omitted are regular in z — w. The transforma- 
tion rule ( IA4I > means that when G (z) is transformed from the 
punctured plane to the cylinder, it is multiplied by a factor of 
z 3 / 2 . Thus the antiperiodicity of the electron/supersymmetry 
generators on the cylinder means that they are periodic on 
the punctured plane. Since we are only interested in this sec- 
tor, the primary fields are those for which the exponent of the 
(z — w) piece in iA5\ is an integer. The lowest non-zero value 
of I7I for which this is true is 7 = l/v3, so we must have 
r = 1/ -\/3. Thus the allowed operators Vj here must have 
7 = n/v3 for integer n. In this language, the supersymme- 
try generators themselves correspond to n — ±3. Thus the 
only primary fields under the extended algebra have n = 
(the identity field), and n = ±1 (which have dimension 1/6). 
The remaining values of n can be obtained by acting with G ± 
on these three primaries, as is clear from the operator product 
expansion \A5l . For example, one obtains V 2 /^ by acting 
with G+ on V_^i^. Physically, this means that all the edge 
modes can be obtained from these three by attaching a hole or 
an electron. There are thus three types of quasiparticles in the 
Laughlin state at ^ = 1/3, so T> 1 / 3 = y3. For general integer 
m = 1/v, the extended symmetry generator is V^^, so by 
the same type of argument, one finds m primary fields. This 
yields 



display here the characters for the free boson. This also will 
provide a check on the above result for T>. 

For a free boson, the highest-weight states under the sym- 
metry algebra including the U(l) current dip are created by 
the primary fields V 1 defined in \A2\ . The highest-weight 
state created by V 1 has dimension 7 2 /2, so the leading term 
in the character must have exponent 7 2 /2 — 1/24. The full 
characters with this symmetry algebra are (see e.g. Ref. 26) 



q P / 4 _q(p+ 2 > / 4 



l=p/V2 



(A7) 



for p integer. The denominator is the Dedekind eta function, 
which is defined as 



q n )- 



Its key property is that when g^ 1 / 24 ^ is expanded in powers 
of q, the coefficient of q m is the number of partitions of m, i.e. 
the number of ways m can be written as the sum of positive 
integers. This often occurs in characters, because descendant 
states are given by acting on primary states with the L_ n s for 
n a positive integer. 

The characters for the symmetry algebra extended to in- 
clude supersymmetry are a sum over the characters X m /</3> 
because the supersymmetry generators V+ ^/g acting on 
ta k e it to V( m ±3)/s/3' as snown m (IA5> . To avoid con- 
fusion with the Ising characters and with the x 7 , we denote 
the character in the identity sector here as x[0], and those for 
the primaries with m = ±1 as ^[±1/3]. One then finds that 



£>„=l/m — 



(A6) 



as noted in the introduction. 

One important subtlety to note is that this is the sector of 
the theory corresponding to antiperiodic boundary conditions 
around the cylinder for fermionic fields such as the supersym- 
metry generator. In the language of supersymmetry, this is 
called the Neveu-Schwarz sector. As we noted previously, 
periodic boundary conditions on fermions in an edge theory 
arise from a vortex in the bulk. Fields local with respect to a 
supersymmetry generator with periodic boundary conditions 
are in the Ramond sector. In an AT = 2 supersymmetric the- 
ory, there must be the same number of primary fields in the 
Ramond sector as the Neveu-Schwarz sector, so D remains 
the same. For the boson with r = l/v3, the fields V n ,^ 
for half-integer n are local with respect to G± with periodic 
boundary conditions. The Ramond primaries are then those 
with n = ±1/2 and n — 3/2, so T> indeed remains V^. Note 
that n = —3/2 is not a primary because it can be obtained by 
acting on V^/ 2 with G-. 

Models with non-abelian statistics need not have d a = 1 
like the free boson does. Thus one needs to work out the fu- 
sion algebra, which is simplest to do by using the modular 
S matrix. Since many interesting models (e.g. all models of 
the fractional quantum Hall effect) involve a free boson, we 



[0] = E x n 



V3 



J(q 3 



71— — OO 
OO 



±1/3] = I J2 (Xn/y/S - XnVS) 



J(q 1/3 ) ~ J(q 3 ) 
277 



where we have defined 



J(p) ^ £ p n2 '\ 



so that for small p, J(p) = 1 + 2p 1/2 + 2p 2 + It is a 

straightforward exercise to generalize these free-boson char- 
acters to the v = 1/m case, where the extended-symmetry 
generators are of dimension m/2. 

To find the modular S matrix, we need to rewrite the char- 
acters as functions of q. Luckily, doing such modular trans- 
formations is a well-understood branch of mathematics. We 
have defined J(p) as a special case of a Jacobi elliptic theta 
function. Then, from standard mathematical texts, one finds 
that 



J(q 1/A ) 

v(q) 



(-Alng/2^) 1 / 2 J( g A ) 
(-lnq/2n)^ 2 r 1 (q) . 



(A8) 
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A little algebra then gives the modular transformations 



X[0](q) = 



x[±l/3]® 



1 



V3r)(q) 



J(g 1/3 ) , 



_V3_ 

Mq) 



2V3r/(g) 



J(q 



We now need to rewrite the right-hand-side in terms of \[0] (q) 
and x[±l/3](g). We can exploit the equality x[l/3](g) = 
x[— l/3](g) to find the modular S matrix which is both sym- 
metric and unitary. We find 



S v 



-.1/3 




From this modular S matrix we can read off the quantum 
dimensions d a = 1 for all three quasiparticles, and T> — 
l/S'o = y3, in agreement with our earlier results. 



3. Moore-Read state 

The edge theory for the v = 5/2 Moore-Read state con- 
sists of a free boson (usually called the charge mode), and a 
Majorana fermion (usually called the neutral mode) 39 . This 
conformal field theory has c = 3/2. The edge electron/hole 
creation operators (i.e. the supersymmetry generators) are 



G d 



ip e 



±iV%<p 



Because the fermion ip has dimension 1/2, the generators 
have dimension 3/2, and it is simple to show that as in the 
case of the f = 1/3 Laughlin state, the extended symmetry is 
J\f = 2 supersymmetry. 

The primaries of the combined field theory must be prod- 
ucts of the primaries /, ip, a of the Ising model with the pri- 
maries Vy for a free boson. The Ising primary fields have 
dimensions 0,1/2,1/16 respectively, so by using the OPE 
(IA3i and the Ising fusion rules (|8}, it is simple to work out 
the powers of z — w in the operator products in the combined 
theory. Demanding that the electron/superconformal genera- 
tors be periodic on the punctured plane gives the primaries of 
the extended symmetry algebra to be 



/, aV 



±1/(2^2)' V ±l/V2- 



(A9) 



All primaries are all invariant under the combined transfor- 
mations a — > —a, ip — > ip and ip — > ip + 2tt\/2\ this is the 
condition that ensures locality. Fields that are odd under this 
symmetry comprise the Ramond sector. 

We can compute the quantum dimensions without the full 
computation of the modular S matrix. The reason is that in 
( IA9> we have decomposed the operators into products of Ising 
fields with bosonic operators, and we know the fusion algebra 
of both. Since the boson and the Ising theory are independent, 
the quantum dimension of the product of fields from the two 
is just the product of the two quantum dimensions. Namely, 



/, ip and Vy have d] = cL 



1, while d a = y/2. The 



total quantum dimension for the Moore-Read state is thus 



V 



MR 



2 + 1 + 1 + 1 = 2V2. 



(A10) 



To check this, we find the characters for the fields in the 
Moore-Read state. These involve both the Ising characters 
and the free-boson characters, but are not just simple prod- 
ucts. The reason is that we want only the states in the an- 
tiperiodic (Neveu-Schwarz) sector, which are invariant under 
ip — > <p + 27r\/2, a — > —a. We want to include descendants 
which are found by acting with the supersymmetry generators 
G , which is invariant under the symmetry. The characters 
for the six primary fields listed in jA9> are then 



aV_ 



±l/2\/2 
V ±l/V2 



XI 



J(q S 



Xip- 



J{q 2 ) ~ J(q 8 ) 



7] 7] 

J(q 1/8 ) - J(q 1/2 ) 
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J(q 8 ) ^ J{q 2 ) - J(q S ) 

Xi> 1- Xi 

7] 7] 

J(q 1/2 ) ~ J(q 2 ) 



(XI + XV ) 



2r, 



Using the Ising modular S matrix and the modular transfor- 
mations in (IA8I one can work out the Moore-Read modular S 
matrix. For the quantum dimensions, we only need the first 
row, which is 



/ 1 11 1 



1 



1 



S " V2V2'2'2'2^'2V2'2V2 



We recover of course the individual quantum dimensions of 
y2 for the two quasiparticles involving the a field, and 1 for 
the other four, as well as the total quantum dimension T>mr = 

1/S§ = 2V2. 

Just to test the formalism a little more, let us compute the 
quantum dimensions for the Moore-Read state for a system 
of bosons at v = 1. The edge theory is very similar, except 
now the extended symmetry is not supersymmetry, but rather 
SU(2). The conformal field theory is still an Ising model and 
a free boson, except now the boson is at a different radius. 
At this radius, the boson can be fermionized into a free Dirac 
fermion. This can be split into two free Majorana fermions 
which, combined with the Majorana fermion from the neu- 
tral sector, form an SU(2) triplet 40 . This model is equivalent, 
via non-Abelian bosonization, to the SU(2)2 Wess-Zumino- 
Witten model, with the extended symmetry being the corre- 
sponding Kac-Moody algebra. The extra SU(2) generators 
are 

J± = V e ±itp 

while J z is the U(l) current dcj). Since the theory is bosonic, 
the symmetry generators are periodic on the cylinder; because 
J± and J z have dimension 1, they are periodic on the punc- 
tured plane as well. The primary fields are therefore 

I, ctVl/2, V 
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The total quantum dimension of this theory is therefore 

V su{2)2 = VI + 2 + 1 = 2. 

Note that a V- 1/2 is not primary because it can be obtained by 
acting with J_ on a V1/2. This is just a fancy way of saying 
that the two form a doublet under SU (2). Likewise, the fields 
ip, e lv , e~ ilp form a triplet under SU{2). 



4. Read-Rezayi 

For a level-/cth bosonic Read-Rezayi stated, the edge 
conformal field theory is the Wess-Zumino-Witten model 
SU(2)k- This has been well studied in the literature, and the 
modular S matrix is given by 41 



S b = 



2 . /7r(a + 1)0+1) 
■ sin 1 



k + 2 V k + 2 
for < a, 6 < fc. The quantum dimensions are therefore 

sin(7r(a + l)/(fc + 2)) 



da sm{w/(k + 2)) 

and the total quantum dimension is 



V 



SU(2) k 



y/k + 2 



V2sin(7r/(fc + 2))' 



(All) 



(A12) 



For the Read-Rezayi fermionic fractional quantum Hall 
states, we cannot instantly read off the modular S matrix from 



the literature. The reason is that we need the characters for 
the extended symmetry in the Neveu-Schwarz sector of the 
J\f = 2 supersymmetric minimal models, which do not seem 
to be in the literature. However, we can compute the quantum 
dimensions by utilizing the fact that the JV = 2 primaries are 
closely related to SU(2) k fields. Namely, the SU{2) k fields 
are labeled by their SU(2) quantum numbers, a spin 1/2 with 
I = 0, 1, ... k and a charge m = — I, —l + 2,...,l. The pri- 
mary fields in the Neveu-Schwarz sector of the superconfor- 
mal field theory are labeled in the same way, and differ only in 
multiplying by a field involving the charge boson, namely V~ 
with 7 = —to/ y/2(k + 2)42^ This changes the field's di- 
mension, but not the fusion rules, because to is a conserved 
quantum number in fusion (in the Hall effect the physical 
charge of the quasiparticle created by the field is me/ (k + 2)). 
Thus a given field's quantum dimension depends only on the 
index I, and all the quantum dimensions of the Read-Rezayi 
fractional quantum Hall states are given by There are 

more primary fields in the supersymmetric case, because in 
the SU (2)fc case, the I + 1 fields with a given I form an SU (2) 
multiplet, and so correspond to only one primary field of the 
S'L/(2)-extended algebra. The total quantum dimension for 
the kth Read-Rezayi state for the fractional quantum Hall ef- 
fect is therefore 



V 



RR 



sin 2 (^ + l)/(fc + 2)) 
sin 2 (7r/(fc + 2)) 



1/2 



2sin(7r/(fc + 2)) 
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